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Abstract 

Using the spinon picture we derive an integral representation for the exact n-spinon 
dynamic structure factor of the spin-1/2 Heisenberg model. 
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1 Introduction 



Since the pionnering work of Bethe |], the Heisenberg model in the anti-ferromagnetic (AF) 
regime has received considerable attention ||2fl- [[19fl . It is also of great experimental interest 
where the dynamic structure factor (DSF) plays a central role in magnetic neutron scattering 
experiments 0. ^TJ. We extend the recent derivation of the 2-spinon DSF of the s — 1/2 AF 
Heisenberg model to the exact n-spinon one [p^] . For an extensive analysis of this 2-spinon 



DSF see Refs. [33, 2"1| . Though more complicated, the result of this paper can be thought 
of as the analogue of that of the s — 1/2 XX model at zero temperature, and which was 
first derived by Niemeijer and later re-derived in a very simple manner by Miiller at al 
nfl. Our work is based on the spinon picture of the eigenspace of the Heisenberg model in 



the AF regime, which was rigorously established in Refs. || |T(|. We use the form factors 
that have been derived in Ref. [|16j from the isotropic limits of those of the XXZ model. 

This letter is organized as follows: first, we briefly review the spinon picture of the 
eigenspace of the s — 1/2 XXZ Hamiltonian. Then, we define one component of the DSF 
for this model. By taking the isotropic limit, we derive an exact integral representation for 
all non-vanishing components of the n-spinon DSF. The summation over all these n-spinon 
contributions leads to the total DSF of the AF Heisenberg model. 



2 The spinon picture in the XXZ model 

The Hamiltonian of the anisotropic (XXZ) Heisenberg model is defined by 

1 oo 

H xxz = -- £ + « + i + A« + i), (1) 

n=— oo 

where A = (q + q~ 1 )/2 is the anisotropy parameter. Here a^' y,z are the usual Pauli matrices 
acting at the n th position of the formal infinite tensor product 

W = ---V® V® V---, (2) 

where V is the two-dimensional representation of U q (sl(2)) quantum group. We consider 
the model in the anti- ferromagnetic regime A < — 1, i.e., — 1 < q < 0. Later we take the 
isotropic limit q — > — 1 in a special manner. This model is symmetric under the quantum 
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group U q (sl(2)), and therefore its eigenspace is identified with the following level U q (sl(2)) 
module P~5[] : 

F = J2V(A i )®V(A ] y. (3) 

Here A$ and V(Aj); z = 0, 1 are level 1 f/ g (s/(2))-highest weights and C/ ? (sZ(2))-highest weight 
modules, respectively, whereas V(Aj)* are dual modules defined from V(Aj) through the 
antipode. The module V(Aj) is identified with the subspace of the following formal semi- 
infinite tensor product of Vs: 

X = ■ ■ ■ V <g> V <8> V, (4) 

and which consists of all linear combinations of spin configurations with fixed boundary 
conditions such that the eigenvalues of a* are (— l) 4+n in the limit n — > — oo. Similarly, 
the module V(Aj)* is associated with the right semi-infinite tensor product of the Vs. The 
spinon picture of the eigenspace of this Hamiltonian is given in terms of vertex operators 
(of type II according to the terminology of Ref. |IJ|) which act as intertwiners of U q (sl(2)) 
modules, and which create the set of eigenstates (i.e., spinons) {|£i, • • •£„ > ei) ... eni i, n > 0}. 
Here % — 0, 1 fixes the boundary conditions, £j denotes a spectral parameter living on the 
unit circle, and 6j = ±1 is twice the ^-component of the spin of a spinon. The actions on T 
of Hxxz arid the translation operator T, which shifts the spin chain by one site, are given 
by 

n 

T\£i,---,£n>i = n^)- 1 ^, •••,£„>!_,, T|0 >,= |0 >!_,, 

i=l 

Hxxz\il,-- -,€n >i = e(£i)|£i, ■ • • , £n >h 



(5) 



where 



(6) 

= i 2f^log^(0 = ^sinh(^)dn(^a). 
Here, e(a) and p(a) are the energy and the momentum of the spinon respectively, am(i) 
and dn(x) are the usual elliptic amplitude and delta functions, with nome — q and complete 
elliptic integrals K and K' , and 

q = — exp(— irK'/K), 
£ = ie iQ , 

(7) 

6 x (y) = (x;x) 00 (y;x) 00 (xy 1 ;x) cx> , 

= n^ =0 (i-^ n )- 



Therefore, a x,v ' z (t, n) at time t and position n are related to a x,y,z (0, 0) at time and position 
through: 

o*»>*(t,n) = e^tH xxz )T- n a x > y > z {W)T n eM-ttH xxz ). (8) 



Moreover, the completeness relation in T reads as [16 



3 An integral representation of the n-spinon DSF of 
the Heisenberg model 

Let us first recall the definition of one of the components of the DSF in the case of the XXZ 
model in the sector % as: 

/oo 
dt e i(wt+km \ < 0|or+(t,m)a-(0,0)|0 >»= ^ Sl; + -(w,k), (10) 

where w and fc are the energy and momentum transfer respectively. It is sufficient to consider 
just the above component because in the isotopic limit, which is our main purpose later, all 
non- vanishing components are equal and proportional to S h ~* (w, k). Due to the conservation 
of the total spin, only pairs of spinons contribute to the DSF. Using the completeness relation, 
the n-spinon contribution is given by 

Sii + -(w, k) = % E me z£ £1 ,..., £n § ■ ■ ■ * A ex P H fc + z*=ip&))) s ( w - £"=i e ^)) 

X i+m < 0|a + (0,0)|f n ,...,fi > en ,...,e 1 ;i+m i;e u ...,e n < £l, ■ ■ ■ , £nW~ (0, 0) |0 > { 
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which, in turn, can be re-written as 



S^iw, k) = 5 E £ll ..., en § ^ - - - § ^ E m6 z exp [2vm{k + Z] =1 P&)\ 
x5(w - E"=i e(Q) (i < 0\a+(0, 0)|f n , ■ ■ ■ , £1 >e n ,...,e 1 ;i i;ei ,„. (6n < £1, • • • , £nk~(0, 0)|0 >, 
+exp (z(fc + E"=i P&)) i-,. < 0|a+(0, 0)|^, . . . , & >,, , : , , i;ei ,..., e „ < 6, • • • , £>"(0, 0)|0 >.. 
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The non- vanishing form factors have been computed in |16[, and satisfy the following rela- 
tions: 

i<0|<7-(0,0)|£ n ,...,£i > e „,..., £i; i = i_i< 0|(7+(0,0)|£n,---,6i >_ e „,..._ ei;1 _i 

= ,< 0|or+(0, 0)| - q£\ . . . , -gC 1 >-e 1 ,...,-e ni <, ( 13 ) 

i;e ll ...,e n <6, • • • , £nk~(0, 0)|0 >, = *< 0|<T _ (0, 0)| - q£ X , -q£ n >- ei ,...,e n ;i ■ 

Their isotropic limits as q — > —1 are performed by first making the following redefinitions: 

c ■ * A 

t = ie™ , 

(14) 

q = -e" e , e -> + , 

with /3, the appropriate spectral parameter for the Heisenberg model, being real. 
Then, one finds the following exact isotropic limits |TO]:f| '0 

U< 0|or+(0,0)|^,..-,ei >e n ,...,e 1 ;i \^---^^ 
2 3 "( 1 -")/ 2 | gei ,..., £n (/3 1 ,...A)] 2 n i < J <y<,J^-(^'-/3 j )l 2 ,o 

7r (7^-14r l +8)/4 r ( 1/ 4 ) n(n-2)| j4 _ (i7r/ /2)|n(n-2) C Osh(/3j ) Pl ' " " P ™ ' 

Iim( 1 _ i < 0|a+(0, 0)|e„, ■ ■ ■ , 6 > en ,.., ei; i-i) = - lim(i< 0|<7+(0, 0)|f n , . . . , 6 >e n ,...,er,i), 
P(Zj) ~> P(&) = cot-^sinh^)), -7T < p(&) < 0, 
e(Q - e{Pj) = = -7rsin(K/3;)), 1 < J < n, 



(15) 



where 



9€-i.,...,fJfix, ■■-,Pn) = ( s J2"=i £ j~ 2 ^ §C t ^Ti \~\i<j<i<n(ai ~ Pj + f ) 

x ni<t<i<n(Pi -«* + ?) n^r(-| + ^)r(-| - ^#) 

x ILeL sinh(c^) FU<*'(a* ~ <*l' + ^i) sinh(c^ - a e ). 



(16) 



Moreover, T(x) is the usual gamma function, and each contour Cg encloses only the poles 
an = j3j — (2m + l)m/2,m > 1 and at = f3j + m /2, with 1 < j < n, and is oriented 
anticlockwise. The set L is defined by 

n 

L = {j, s.t. €j = +1, and ^e* = -2, 1 < j < ra}. (17) 

i=l 



1 We do not find the same overall coefficient in this limit as that of Ref. [jl6| 

2 We thank Karbach and Miiller for confirming to us the relation | J 4_(«7r/2)| 2 1 ^4+ (*7r/2) | 2 = 1/2, which is 



used here. See [[16) for the definition of A+(x) 
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Restricting to the first Brillouin zone, that is, < k < 2tt, and performing the two delta 
integrations we find that S^ + ~(w, k — tt) is independent of i and reads for n > 2 



S^-(w,k-n) = C n &dp 3 ...J%dp n Efr&)Q(W u -W)Q(W-Wi) 

/(/3(pi),/3(p 2 ),/3(P3),-,/3(Pn))X; eii £n \ge 1: ...,, n (0(pi),P(p2),f3(p 3 ),...,f3(p n ))\ 2 

Here, 6 is the Heaviside step function and 

f(P(Pl),P(P2),P( P 3), ■ ■ .,P(Pn)) = Ul< 3<r <n \A-(PM ~ P(Pj))\^ 

„ 2(- 3n2 + 2n + 4 )/ 2 7r 7n ( 2_n )/ 4 
<~<n = „!r(l/4)"(™- 2 )|A_(i7r/2)|' l (™- 2 ) ' 

W = w + ttJ:] =3 unfa), 
K = k + E] =3 Pj, 
W u = 2rr\sm(K/2)\, 
Wi = 7r| sin K\. 

Moreover for fixed W and K, the sum S(pi, P - 2 ) is over ah pairs (pi,P2), solutions to the 
energy-momentum conservation laws: 



(19) 



W = -7r(sin(pi) + sin(p 2 )), 
K = -p 1 -p 2 . 



(20) 



The total n-spinon contribution is obtained through 

S+-(w, k - tt) = £ fc — 7r) = 2S£+->, fc - 7T), (21) 

i=0 

from which we derive all the non-vanishing components of the DSF as: 

SW{w,k-n) = 2S+-(w,k-ir), fx = x,y,z. (22) 

The 2-spinon case, i.e., n = 2, is particularly interesting because all contour integrals and the 
complicated g function in (|18|) drop out. Furthermore, because just the four pairs (pi, p^), 
(p2,Pi), (—Pi — — pi — tt) and (— p~ 2 — tt, —pi — tt) contribute, and moreover equally, to the 



sum Z)(pi,p a )j £2 simplifies drastically to |22 



S2~{w k — ir) = QC. z ( w u-- w ) @ ( w - w i)\ A -(f 3 (pi-)-f 3 (P2)\ 

v^"™ 2 ' (23) 



e(w„-w)e(w-w;)|A-(/3(pi)- / 3(p- 2 )| 2 



where now 

w = — 7r(sin(pi) + sin(p 2 )), 

^ ^ (24) 
w u = 27rsin(/c/2), 

u>z = 7r I sin | . 

Let us mention that S2~(w, k — n) has been expressed in Ref. |23| just in terms of the energy 
and momentum transfer through the following relation: 



\ 



w 2 u - wf 
w 2 — w? 



(25) 



It is also interesting to note that 5*^ (w, k — n) vanishes at the upper boundary w = w u 
since ^4_(0) = 0. This is consistent with both facts that at the latter boundary pi = p 2 
and that the spinons obey the Pauli exclusion principle common to fermions. More detailed 
analysis of the singular behaviour of S2~(w, k — n) at both boundaries w = w u and w = wi 
can be found in Ref. [ 23|| . Finally, let us mention that it would be particularly interesting to 
evaluate the contour integrals that appear in the g function of relation (|18[), and especially 



for n = 4. The latter case is also useful in optical experiments p6|. 
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